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Abstract 



Resonant formation of d/it molecules in collisions of muonic tritium (/it) 

on D2 was investigated using a beam of /it atoms, demonstrating a new direct 

approach in muon catalyzed fusion studies. Strong epithermal resonances in 

d/it formation were directly revealed for the first time. Prom the time-of- 

flight analysis of 2036 ± 116 dt fusion events, a formation rate consistent with 

0.73 lb (0.16),„eas i {^■O^)model times the theoretical prediction was obtained. 

For the largest peak at a resonance energy of 0.423±0.037 eV, this corresponds 

to a rate of (7.1 ± 1.8) x 10^ s~^, more than an order of magnitude larger than 

those at low energies. 
36.10.Dr, 21.45.+V, 25.60.Pj 
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Typeset using REVT^ 



Reactions of muonic hydrogen atoms and molecules present a sensitive testing ground 
for few-body theories involving strong, weak, and electromagnetic interactions. Among 
such reactions, muon catalyzed fusion (yuCF) in a deuterium-tritium mixture has attracted 
particular interest, where one muon can catalyze more than 100 nuclear fusions between 
deuteron and triton {d + t a + n) ^j^. A key step in high-yield /xCF is resonant 
formation ^ of the dfit molecule 

/it + [{dfit)ndee]Ku, (1) 

where the collision energy, and the energy released upon formation of dfit in the loosely 
bound state {Jv = 11), are absorbed in the ro- vibrational {Ku) excitation of the molecular 
complex [{dfit)dee], a hydrogen-like molecule with {dfit)~^ playing the role of one of the 
nuclei. Because of its limiting rate Xd^t, the process has been considered one of the major 
bottlenecks in achieving high efficiency in fiCF. However, theoretical calculations predict 
strongly enhanced resonances for reaction ([1|) at fit kinetic energies of order 1 eV (epithermal 
resonances or ER), yet the experimental information is scarce thus far. It is the purpose of 
this letter to report a direct confirmation of the predicted ER and a determination of their 
energies, enabled for the first time with the atomic beam approach [Q. 

In conventional /iCF experiments with homogeneous mixtures of hydrogen isotopes 
the extraction of A^^t and its energy dependence relies on a kinetic model describing a 
complex chain of reactions, which includes processes that are not well understood, and 
assumptions that have been recently challenged (see Ref. 0). The collision energy for 
reaction (1) in the equilibrium states is given by the target thermal energy, but difficulties 
in tritium handling have so far prevented realization of high temperature targets capable of 
thermally accessing ER. Although transient phenomenagave first evidence and subsequent 
insight into the epithermal effects in conventional targets, where a small fraction of dfit 
formation may occur before the thermalization, no quantitative measurement of Xd^t for ER 
or its energy E^^g licis been obtained due primarily to unknown fit initial conditions. 

The atomic beam method described here provides direct access to ER because of the 
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available fit beam energy (0.1-10 eV). A single dfit formation can be studied, on an event- 
by-event basis, isolated in time and space from other processes. With the fit time of flight 
(TOF) between the separated layers providing a measure of the collision energy, Eres can 
be determined in a direct manner. Furthermore, manipulation of heterogeneous multi-layers 
allows us to study and control, individually, important processes such as fj, transfer, fit 
emission, and moderation. 

The experiment was performed at the M20B channel at TRIUMF. Details of the mea- 
surement and analysis may be found in Ref. |10|. The apparatus [p!li|l2| is illustrated in 



Fig. 1. Target layers (Fig. 2(a)) were prepared by rapidly freezing a hydrogen isotope or 
mixture onto the gold foils, held at 3.5 K in an ultra-high vacuum of order 10~^ torr or 
better. The layer thickness and uniformity were characterized off-line via energy loss of a 
particles , and the effective thickness was then derived using a separate measurement |13 



of the muon beam profile (about 30 mm diameter). 

A beam of 5 x 10^ fi^ s^^ of momentum p=27 MeV/c and momentum spread Ap/p=5.5% 
(FWHM), defined by a 250 fim scintillator of diameter 48 mm (Tl), entered the target 
vacuum. It was degraded mainly by a 51 fim gold foil upstream (US) from a TOF region. 
A fraction Sp of the fi stopped in a tritium-doped hydrogen layer (an emission layer) of 
3.43±0.18 mg-cm~^ frozen to the US foil; most of these initially formed fip. The transfer 



fip — i> fit took place in a time of typically 100 ns [Tj], creating fit with recoil energy of 45 
eV due to the reduced mass difference. Because of the Ramsauer-Townsend effect in fit + p 
scattering, fit atoms are emitted with energies near 10 eV into an adjacent layer |13|. A 



D2 moderation layer (96.0±5.0 /ig-cm~^) efficiently reduced the fit beam energy to ~1 eV 
via elastic scattering in order to better match ER energies. The fit, after a flight time of a 
few fis, reached the D2 reaction layer downstream (DS) from the TOF region (17.9±0.5 mm 
from the US layers) and formed [{dfit)dee] from which fusion can occur to produce a and n. 
The time between the fi entrance signal and the detection of a fusion a (which we call fusion 
time) is dominated by fit TOF and provides information on molecular formation energy, as 
long as the energy loss (AE^t) of fit, due to elastic scattering before the formation in the DS 



D2, is small. A thin DS layer of 21.2±1.4 /ig-cm~^ was chosen to minimize ^E^t so as not 
to obscure the time-energy correlation. The detailed theoretical description of our method 
is given in Ref. ||15|| . 

Two series of data, with emission layer tritium concentrations q = 0.1% (Run 1) and 
Q = 0.2% (Run 2), were analyzed separately. Deuterium concentrations in H2 were typically 
less than 2 ppm. Run 1 had a 1.77±0.12 mg-cm^^ H2 substrate beneath the D2 DS reaction 
layer, while for Run 2 the D2 reaction layer was deposited directly on the DS gold foil. 
The data were normalized to the corrected number of incident muons, A*"^, which takes 
into account the data acquisition dead time (~20%) and pile-up of incident muons in the 
10 jjis gate (~5%). The signal and background in the Si detectors for Run 1 are shown 
in Fig. 2(b) with a delayed time cut selecting DS events, giving a S/B ratio of about 2:1. 
The background, mainly due to /x decay and capture related processes (few fusion events 
from the US D2 moderator pass the delayed time cut because there is no TOF delay), 
can be accurately determined with the beam method from runs without the DS reaction 
layer, in which only DS fusion is turned off while other processes are not affected. Potential 
background from protons from /i-induced dd fusion, arising for example from a recycled 
following dt fusion, were estimated to give a 2.4 ± 0.9% correction in DS fusion by looking 
at a-p correlations in the Si detectors, using the much larger data sample from fusion in 
the US moderation layer. A small contamination (^ 5 ppm) of nitrogen in the target layer 
was estimated to reduce the fusion yield by about 2%. Any further residual effects {e.g., 
due to inaccuracy in the scaler, or time zero shifts) were carefully investigated and were 
conservatively reflected in the final errors. With totals for of 6.02 x 10*^ and 2.82 x 10*^ for 
the data and background runs, respectively, we have observed 2036 ± 116 DS fusion events, 
for Runs 1 and 2 combined. 

An absolute measurement of the fusion yield required determination of several factors 
including stopping fraction Sp, Si detector solid angle VLsi, and energy and time cut accep- 
tances, eE and e-r- Our value of Sp = 0.299 ± 0.015 was based on fits of the decay electron 
time spectra recorded by electron and neutron counters (in the charged mode), where n 
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stopped in H2 exhibited a characteristic ~2 /is hfetime. The absolute efficiency of the coun- 
ters was determined from the detection of delayed electrons following the observation of 
fusion Sf then could be derived model-independently from the absolute amplitude at 
time zero of the 2 /is decay electron component, normalized to A*"^. An independent and 
consistent estimate of Sp was obtained from a GEANT |TB] beam and decay simulation, 



which reproduced our measured range curve [|I7 



The energy cut efficiency and possible existence of tails in the a energy distribution in the 
Si spectrum (Fig. 2(b)) were investigated by changing the cut width, as well as by simulating 
the energy loss of a particles in the DS D2 layer giving ep = 0.978±0.064 for the applied 
cut of 3.1 < Ea < 3.7 MeV. ilsi was determined to be 2 x (2.46 ± 0.10)% from similar 
calculations, in which the spatial distribution of fusion events was constrained by electron 



imaging of muon decay. Effects of possible beam shifts [|T9[ as well as small geometrical 
inaccuracy are reflected in the errors. After combining all the factors, including a small time 
cut correction, the absolute normalization was determined to the relative precisions of 9.5% 
and 7.7% for Runs 1 and 2, which together with statistical uncertainties gave the DS fusion 
yield per stopped of (5.04 ± 0.50) x lO^V/x and (4.27 ± 0.54) x lO^V^u for Runs 1 and 2, 
respectively. 

The results were used to test theoretical predictions as follows. The formation rate X^J'^^ 
from Ref. Q, corrected for the Doppler broadening due to the target molecule motion at 
3 K (without many-body effects), was used as a standard input, with low energy rates 
further modified to reflect possible sub-threshold effects 0. The experimental TOP fusion 
distributions were ffited with a series of Monte Carlo (MC) simulated spectra for which the 
input Xd^t and Eres were varied as Xd^t = 'SxX^Jp.t Eres = SeEI^^. The measured absolute 
normalization and its estimated uncertainty were used to constrain the fit. The best values 
for S\ and Se, which were taken to be energy- and rate- independent, respectively, were 
extracted from the distribution. The fusion probability Wj from Ref. was kept fixed 
during the ffi. A detailed simulation code was developed for the analysis |2^. Input cross 
sections for muonic processes were based on theoretical values in Refs. |pTlj2^ , but our multi- 
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layered target allowed independent tests of important reactions including fip — > fit transfer 



and pfip formation [111], and fit + p and fit + d collisions |T0|,|T3|. In addition, independent 



simulation codes p!5|j23|] were used to check the consistency of some of the key processes. 
We stress the importance of explicitly treating the resonant scattering in the calculations; 
it can be shown that the use of the renormalized effective rate A^^^ ~ WfXd^t as in 
Ref. [|l5l would result in the overestimate of the fusion yield by as much as a factor of two, 



a fact which resolves an earlier discrepancy reported in Ref. The effect of the energy 
distribution of resonantly scattered fit was investigated within the kinematically possible 



range. Solid state effects in thermalization and formation [^, which are important mostly 
at lowest energies, were not explicitly included in the simulations, but their influence at ER 
energies is expected to be small. The effects were estimated by appropriately modifying the 
cross sections. Table | summarizes our evaluation of systematic uncertainties in the MC 
modelling. 

Figure ^ shows a comparison of the calculated MC spectrum with the data for Run 
1. Also plotted are the simulated contributions from the time-energy correlated events 
{i.e. with small AE^t), exhibiting more clearly the resonance structures due to different 
u states. From the fits, we obtained Sx = 0.88 ± 0.11 (Run 1) and 0.55 ± 0.12 (Run 2) 
with xVdof (dof = 50) of 0.96 and 1.27 respectively, and Se = 0.928 ± 0.040 (Run 1) and 
0.994 ± 0.087 (Run 2) with similar fit quality. The fit uncertainties reflect both statistical 
and normalization errors. An apparent deviation at the 2a level in Sx between Runs 1 
and 2 may indicate some unaccounted-for systematic uncertainties (due possibly to slightly 
different experimental conditions as described above), hence the measurement error was 
increased accordingly so that two data set give of 1 . When the known MC modelling 



error (Table |) is included, the combined final result is Sx = 0.73 ± (0.16)m,eas ± (0.09)mode«- 
As for the resonance energy measurements, the weighted average of two runs gave Se = 
0.940 ± (0.036)meas ± (0.080)mode«, whcre the modelling error is mainly due to uncertainties 
in the fit + d elastic scattering process and the fit TOF drift distance. 

These scaling results correspond, in the theory [Q, to a peak formation rate for the 



main resonance (z/=3) of (7.1 ± 1.8) x 10*^ s"^ at Eres = 0.423 ± 0.037 eV for /i^^=^ and 
(6.5 ± 1.6) X 10^ s"^ at 0.508 ± 0.047 eV for fit^^^, where F refers to the /it hyperfine state 
(the rates are normahzed to hquid hydrogen density, 4.25 x 10^^ cm~'^). The obtained rates 
are more than an order of magnitude larger than those at low energies experimentally 
demonstrating the prospect for high cycling /zCF in a high temperature target of several 
thousand degrees. 

Our TOF sensitivity for the energy dependence allows us to clearly reject, for example, a 
constant Xd/it, establishing the existence of resonant structure in epithermal dfit formation. 
If one assumes the energy level spectrum of the [{dfit)dee] molecule, our results for Eres 
imply sensitivity to the binding energy of the {dfit)ii state with an accuracy of the order of 
the vacuum polarization corrections. 

We have recently collected data for resonant dfit formation in fit + HD collisions, for 
which the predicted resonances are even stronger The results of that measurement will 
be reported in a future publication [^. 
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TABLES 

TABLE L Estimated effects on rate scaling parameter S\ by the systematic uncertainties in 



the MC modelhng. 

MC error source ASx/Sx % 

H beam size 1.2 

Nonuniform /x stopping (GEANT) 1.8 

fjLt TOE drift distance 2.6 

HP — > jit transfer, pjip formation 5.6 

Ht + d, ixt + t scattering, layer thickness 5.7 

jJLt + p RT minimum energy 1.1 

Resonance Doppler widths in solid 8.0 

Solid state and low energy processes 5.0 
(subthreshold resonances, slow thermalization, 
jit energy after resonant scattering) 

Total (in quadraturo) 12.9 
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FIG. 1. Top view of the apparatus. Si detectors were placed in vacuum viewing target layers 
without a window, enabling the high resolution detection of fusion a particles (Fig. 2(b)). The Ge 
detector monitored target impurities via muonic X-rays, while plastic scintillators (Enl, En2, Ege) 
were used both to detect muon decay electrons and to veto charged particles for Ge and neutron 
(Nl, N2) detectors. 
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FIG. 2. (a) Schematic of the target consisting of emission, moderation, and reaction layers (not 
to scale), (b) Si energy spectra showing signal (histogram) and background (hne) with a time cut 
of t > 1.5 /is selecting DS events delayed by [it TOF. 
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FIG. 3. Time-of-flight fusion spectrum (error bars) and simulation spectrum (histogram). Also 
plotted are simulated contributions from different resonance peaks given by time-energy correlated 
events with AE^t < 0.15 eV. Note that angular dispersion of the nt beam also contributes to the 
peak widths. 
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